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The structure of a finite particle cluster is typically determined by total energy minimization. Here 
we consider the case where a cluster of soft sphere dipoles becomes active, i.e. when the individual 
particles exhibit an additional self-propulsion along their dipole moments. We numerically solve 
the overdamped equations of motion for soft-sphere dipoles in a solvent. Starting from an initial 
metastable dipolar cluster, the self-propulsion generates a complex cluster dynamics. The final 
cluster state has in general a structure widely different to the initial one, the details depend on 
the model parameters and on the protocol of how the self-propulsion is turned on. The center-of- 
mass of the cluster moves on a helical path, the details of which are governed by the initial cluster 
magnetization. An instantaneous switch to a high self-propulsion leads to fission of the cluster. 
However, fission does not occur if the self-propulsion is increased slowly to high strengths. Our 
predictions can be verified through experiments with self-phoretic colloidal Janus-particles and for 
macroscopic self-propelled dipoles in a highly viscous solvent. 

PACS numbers: 82.70.Dd, 61.46.Be 


I. INTRODUCTION 

Clusters comprising N individual particles occur in 
widely different areas of physics ranging from the atomic 
world [l| to nanoparticles [2[, colloids and to 

macroscopic granulates [T2|. In the simplest case, the 
particles interact via a pairwise potential, such as a 
Lennard-Jones potential |1^ or a hard-sphere-dipole in¬ 
teraction |l3-[^ , and the equilibrium groundstate struc¬ 
ture of the cluster is obtained by minimization of the to¬ 
tal potential energy. Even for small N the structure can 
be nontrivial and differs substantially from a cutout of a 
simple crystal. Here, we consider clusters where the con¬ 
stituents are active or self-propelled particles which are 
swimming in a viscous solvent at low Reynolds number. 
Such active particles or microswimmers can be artificially 
realized as colloidal Janus particles exposed to a thermal 
gradient [ 2 ^ close to a solvent phase transition a or a 
chemical reactant catalyzed at one part of the Janus par¬ 
ticle [ 2514 ^. This brings the particle into motion along 
its symmetry axis, thereby creating a nonequilibrium sit¬ 
uation where the minimization principle of the potential 
energy breaks down. An “active cluster” composed of ag¬ 
gregated self-propelled particles will therefore exhibit a 
characteristic structure, and a characteristic motion. In 
particular, as the self-propulsion force grows stronger, the 
cluster can break revealing an activity-induced fission. 
The occurrence of clusters within active suspensions has 
been studied frequently in the last years through experi- 
ments theory (sl, [39|, and simulations jiol-li^ 

considering purely repulsive as well as slightly attrac¬ 
tive particle interactions [ij. However, controlled fis¬ 
sion of such clusters has not been a focus of research. 
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with first forrays in this field made only recently by Soto 
and Golestanian ISlii- 

We describe these phenomena by a simple model com¬ 
bining Weeks-Chandler-Andersen and dipole pair poten¬ 
tials with an intrinsic effective self-propulsion force [13, 
im along the dipole moment [49| . Our model differs from 
the one employed by Soto and Golestanian [ill, EH , which 
focuses on clusters of self-propelled particles induced by 
chemotaxis and characterized by their violation of the 
actio-reactio principle. 

In order to keep the model as simple as possible, we 
consider in this paper only a small number N of dipoles 
for the trivial cases A^ = l,2uptoA^ = 5. The particles 
are soft spheres with an embedded dipole moment that is 
oriented along the direction of self-propulsion. The sit¬ 
uation with = 4 and N = 5 dipoles already reveals 
a quite complex dynamical behaviour. Our initial con¬ 
figuration is either a linear or ring-like structure which 
represents the groundstate of the system [S E El 
or a metastable structure such as various compact and 
branched clusters, in particular a Y-junction of Y = 4 
dipoles SEl- We then introduce the self-propulsion 
and follow the dynamics of the particles. When doing 
so, we distinguish between two different protocols: an 
instantaneous switching to a given velocity and an “adi¬ 
abatic” slow increase of the self-propulsion. As a result, 
we find that the cluster’s center-of-mass generally moves 
along a helical trajectory. The details of the helix are gov¬ 
erned by the magnetization of the initial cluster. More¬ 
over, we find a plethora of final cluster states, which de¬ 
pend on the model parameters and the enacted protocol. 
For strong self-propulsion and a nonlinear initial state 
with nonvanishing dipole moment, fission of the cluster 
occurs for instantaneous switching. However, interest¬ 
ingly, for the same parameters, there is no fission in the 
adiabatic switching case. 

Our results can be verified for colloidal Janus parti- 
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cles which a strong dipole moment along their symmetry 
axis. At vanishing self-propulsion these particles have 
been considered quite a lot in the context of ferrogels 
and ferrofluids . Dipolar self-propelled parti¬ 

cles have been prepared recently by Baraban and cowork¬ 
ers [ 53 , an d aggregation into clusters has indeed been 
observed 

Our paper is organized as follows: In Sec-Illwe specify 
our model for active dipolar particles. The different types 
of motion of clusters with and without vanishing total 
dipole moment are presented in Sec. IIIIl In Sec. [IVl we 
study the fission induced by activity in detail and we 
conclude in Sec. m 


II. MODEL 


We consider N spherical dipoles in three spatial 
dimensions. The position of the ith dipole will 
be denoted as = [xi^yi^Zi] and the dipole mo¬ 

ment nii = mui is directed along the unit vector 
= (sin6>i cos^i, sin6>i sin^i, cos6>i). The total pair¬ 
wise interaction potential Uij is the sum of a Weeks- 
Chandler-Andersen which describes a repulsive 

soft core [56| and a point dipole potential 
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where Yij = — rj is the position of particle j relative to 
particle i and Vij the respective distance, see Fig.[TJ Here, 
the cut-off length is Tc = 2^/^cr. We introduce the self¬ 
propulsion by means of an effective internal force = 
Fui dil, [53 which is directed along all dipole moments 
m^, leading to a constant self-propulsion velocity = 
vvii for an individual single particle. The velocity is given 
by = F//t, and ft denotes the translational Stokes 
friction coefficient. As units of energy and length we 
choose the parameters e and a from the WCA-potential. 
Time is measured in units of r = /tcr^/e. 

The motion of microswimmers is restricted to the low 
Reynolds number regime and the corresponding over¬ 
damped equations of motion for the positions and ori¬ 
entations Ui are 


ft ■ dtVi{t) = -VrV + ftviliit), (3) 

fr • dtUi{t) = -Ti X (4) 

where Tj = Ui(i) x Vu47 is the torque on particle i and 
Uij the total interaction potential. The 
rotational friction coefficient fr of the spherical particles 
is assumed to be linked to ft via the equilibrium relation 
fr = /tcr^/3. We neglect any thermal fluctuations and 



FIG. 1. Sketch of a pair of dipoles in three spatial dimen¬ 
sions with center-of-mass distance Vij. The self-propulsion 
and the dipole moment are directed along the unit vector 
Ui = (sin Oi cos <fi , sin Oi sin (j)i , cos Oi ). 


any solvent-induced hydrodynamic interactions between 
the particles in our model. 

Simulations are performed in three spatial dimensions 
using a simple Euler integration algorithm with time 
steps At = 10for simulation times t = 400r, which 
are sufficiently large enough to allow the active dipole 
clusters to achieve steady state structures. The dipole 
strength is varied in the range 0 < m? / (ecr^) < 6. 

Starting configurations are gathered by an energy min¬ 
imization for the respective parameter sets for passive 
dipoles, V = ^. We consider two protocols of how the 
activity is applied to the dipoles. The self-propulsion is 
either instantaneously increased for the starting configu¬ 
rations, or we slowly increase the activity stepwise by a 
velocity increase of /S.v = 0.1 cr/r subsequent after a long 
waiting time of t = 400r. This corresponds to slow or 
’adiabatic’ switching. 


III. ACTIVE DIPOLE CLUSTERS: SIMPLE 
CASES AND CENTER-OF-MASS 
TRAJECTORIES 


A. The special cases N — 1,2,3 

Let us now discuss the simple cases of very few particles 
N = 1,2,3. First of all, a single particle N = 1 will 
trivially move with velocity v on a straight line along its 
orientation u. Next, the ground state of two dipoles {N = 
2) is a head-tail-configuration where the two magnetic 
moments possess the same orientational direction and are 
aligned along the separation vector of the two spheres. 
This is shown in Fig. He). Putting a drive along the 
magnetic moments will lead to joint motion of this linear 
cluster with the same velocity v as that for the individual 
spheres. This behaviour does not depend on the protocol 
of the drive. 

Third, a magnetic triplet {N = 3) will form a 
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metastable ring-like structure where the sphere centers 
are on an equilateral triangle and the orientations have 
relative angle differences of 27r/3, see Fig. [2fb). How¬ 
ever, the groundstate is a linear chain im. This ring-like 
cluster will rotate around its non-moving center-of-mass 
when a self-propulsion is turned on. The distance be¬ 
tween the spheres and the angular velocity increase with 
increasing drive. 

B. Center-of-mass trajectories 

Next, we consider the center-of-mass coordinate of 
dipole clusters as defined by 

1 ^ 

= (5) 

i=l 

Due to the reciprocal interactions between the individ¬ 
ual dipoles, the motion of the center-of-mass Rc can be 
determined by the orientations Ui of the active particles, 

AT f ^ 

^ i=l i=l 

such that the center-of-mass velocity Vc is proportional 
to the total magnetic moment 
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The motion of the center-of-mass can be easily classified 
as summarized in Fig. O The simplest situation is an 
initial ring cluster as shown for A/" = 3,4, 5 in Fig. Efb) 
whose center-of-mass is in the ring center. This ring has 
a vanishing total dipole moment, M = 0. Based on a 
simple analysis of the equation of motions and on sym¬ 
metry arguments, an initial ring cluster will just rotate 
around its center-of-mass such that the center-of-mass 
is non-moving. This is valid for any for any self¬ 
propulsion strength v and for any protocol (similar to the 
special case N = 3 discussed above) [58|. Note, however, 
that the radius of the rotating ring and the correspond¬ 
ing angular velocity do increase with increasing drive v. 
The result of a non-moving center-of-mass can also be 
obtained for other initial clusters with vanishing initial 
M. An example for a compact three-dimensional cluster 
different from a ring is shown in Fig. EJb) for N = 4. 
This tetrahedral cluster is mechanically stable (but ener¬ 
getically metastable) in equilibrium {v = 0) and has no 
spontaneous magnetization (M = 0). 

Obviously, linear chains of A" > 1 dipoles move on 
a straight line and do not change their relative shape, 
see Fig. [2fc), since all dipole moments and consequently 
the self-propulsion velocities are directed along the main 
symmetry axis. There are more complex structures like 
a Y-junction with A = 4 dipoles, see again Fig. Efc), 
whose center-of-mass moves on a straight line. The rel¬ 
ative positions of these cluster particles, however, can 
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FIG. 2. (a) Trajectories of the center-of-mass for different 

initial active dipole clusters composed of A = 1,..., 5 dipoles, 
(b) Ring-like clusters perform a rotation around the center-of- 
mass such that the center-of-mass is non-moving, (c) Chains, 
as well as a Y-junction, move along a straight line and (d) 
a compact cluster of A = 5 dipoles propagates on a helix. 
Please note that the relative positions of the cluster particles 
change during the cluster propagation for an initial Y-junction 
and the compact cluster shown in (d). For all cluster confor¬ 
mations the respective total potential energy Uje is given for 
m^/(ecr^) = 1. 


change. The details of the relative motion will be con¬ 
sidered in the next section. Finally, a moving cluster of 
constant shape generally performs a helical motion [5^ . 
The helical motion is generated by a non-vanishing mag¬ 
netization M 7 ^ 0 which provides the translational force 
and a non-vanishing torque on the cluster. A typical ex¬ 
ample is an initially compact cluster with non-vanishing 
magnetization as shown for A = 5 particles in Fig. [2fd). 
This compact cluster has the dipoles arranged at the cor¬ 
ners of a pyramid. 


IV. REORGANIZATION AND FISSION OF 
ACTIVE DIPOLE CLUSTERS 

We now turn to clusters composed of four (A = 4) 
and five (A = 5) particles. These cases reveal nontriv¬ 
ial and interesting dynamical behaviour caused by self¬ 
propulsion. In particular, some structures will reorga¬ 
nize, reassemble and split due to the self-propulsion and 
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this depends explicitly on the protocol applied to turn 
on the self-propulsion. To be specific we first consider 
an initial metastable Y-junction for Y = 4 particles and 
then study a compact metastable initial cluster for Y = 5 
particles. 


A. Y-junction (Y = 4) 

An initial metastable Y-junction built up by Y = 4 
dipoles is shown in Fig. me). We study its dynamical 
behaviour for varying dipole strength, self-propulsion ve¬ 
locity V and switching protocol. In this case the motion 
is two-dimensional in the plane set by the initial config¬ 
uration. 

Figure Ha) shows the cluster state diagram after a 
long simulation time, t = 400r, in the two-dimensional 
parameter space of dipole strength and an instanta¬ 
neously applied self-propulsion velocity v. As a result, 
seven different final configurations can be discriminated 
by suitable order parameters, as explained in detail be¬ 
low. Clearly, for small self-propulsion the cluster retains 
its shape as it is metastable but exhibits a drift along its 
total magnetization M. This is indicated by the (red) 
filled circle. Obviously one finds this situation when the 
self-propulsion velocity is small. For very large drives 
on the other hand, fission of the cluster shows up 
(black plus symbols) where the separating dipoles fly 
away pointing forward. In between these two extremes, 
there are five other states: i) a Y-junction where the 
orientations of the ’head dipoles’ now point backwards 
(brown filled square), ii) a kite cluster where these dipoles 
shift back towards the center of the whole cluster (yel¬ 
low filled diamond), iii) a configuration which we denote 
as anti-Y-j unction since the former ’head particles’ are 
located in the back (green filled pentagon), iv) a state 
where two clusters occur, the backbone and two paral¬ 
lel dipoles moving in the opposite direction (blue non¬ 
filled squares) and finally another fission state where the 
separating dipoles burst away pointing backwards (black 
cross). The corresponding trajectories of the dipoles dur¬ 
ing rearrangement of the clusters can be seen in Fig. me), 
in the center-of-mass frame of reference [58l |. 

It is instructive to compare the scaling behaviour of 
our findings with that of hard sphere dipoles. Since 
hard spheres only bear a length scale (such as the di¬ 
ameter cr) and no energy scale, the only ratio which mat¬ 
ters, is that of the dipole-dipole interaction force, propor¬ 
tional to relative to the internal driving force F. 

Therefore, in order to test this scaling, we have plotted 
in Fig. [3jb) the same diagram with the scaled dipolar 
strength vn ?If hard sphere scaling would hold, 
the transition lines should be horizontal in this plot. De¬ 
viations from this scaling then need to be attributed to 
the softness of the WCA-potential. Indeed, as is revealed 
in Fig. [31(b), the transition lines between different states 
are almost horizontal. The largest deviation occurs for 
the transitions where the kite-like structure (yellow filled 


diamond) and the backward Y-j unctions (brown filled 
square) are involved. 

Next, we change the protocol to adiabatic switching, 
see Fig. |3Kd). By slowly increasing the self-propulsion 
velocity, the initial Y-junction deforms into a kite cluster. 
Further increase of v leads to temporal detaching from 
the backbone and formation of an anti-Y-j unction, but 
fission does not occur. This can intuitively be attributed 
to the fact that the attraction of the dipole cluster has a 
larger impact on the cluster connectivity if the particles 
are perturbed smoothly over time. 

We now discuss the details of the different structures 
with the help of suitable order parameters. Let us first in¬ 
troduce the positional angle a of one of the ’head dipoles’ 
position relative to the backbone of the Y-junction, see 
Fig.HKa). We further define an angle f of the orientation 
of the ’head dipole’ relative to the backbone, see again 
Fig. IHa). In general, the angles a and (j) are different. 
A third order parameter is provided by the distance d 
of the head-dipole to the next backbone particle. We 
define connectivity of the cluster by the simple distance 
criterion: if d < Vc the cluster is connected, else it is 
split. Figure SKb) shows the final angle a for fixed dipole 
strength and varied self-propulsion velocity v in the case 
of the instantaneous switching protocol. In line with the 
state diagram discussed previously, for increasing v, the 
angle a broadens and the initial Y-junction {a < 90°) 
becomes a kite cluster (90° < a < 120°), which reorga¬ 
nizes into an anti-Y-j unction (120° < a < 180°). Here, 
the nearest neighbour is now the dipole in the rear of 
the backbone, before finally fission occurs. In the limit 
of small dipolar interactions, an arbitrary self-propulsion 
leads to fission, see again Fig. |3l(a). The test of hard- 
sphere scaling, i.e. a scaled plot versus m^/(Fcr^) [see 
Fig. |3Kb)], reveals a Master curve up to m^/(ecr^) < 3. 
This shows that the hard sphere limit (e ^ oo) is reached 
quickly and that the state diagram essentially depends on 
the parameter Fa^jw? = jw? for m‘^ j< 3, 

which is consistent with Fig. |3l 

Likewise, in Fig.|l|(c) the angle f is shown, which pro¬ 
vides a better resolution of structural details. The ini¬ 
tial passive Y-junction has an angle fo = 53.4°. Self¬ 
propulsion first induces a broadening of f in the Y- 
junction in the range (0o < 0 < 90°) (red filled cir¬ 
cles). Then the orientation of the outermost particles 
flips to the opposite direction (90° < (j) < 180°) (brown 
filled squares) until the kite cluster occurs, where the 
head dipoles turn until they are reversed relative to the 
backbone magnetization, (j) « 180°. This conformation 
only shows up for small dipolar strengths m?. For even 
larger the Y-junction temporarily splits into a remain¬ 
ing backbone and two single dipoles. These dipoles will 
be attracted again to the rear of the backbone, and reat¬ 
tach, being orientated again in the direction of motion 
(270° < (j) < 360°). This reassembled structure is de¬ 
noted as an anti-Y-j unction. If v is increased further, 
these detached dipoles may still be attracted by the back¬ 
bone. However they may just collide in its wake and align 
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FIG. 3. (a) Emerging state diagram for an initial Y-junction for varied reduced dipole strength m^/(ecr^) and self-propulsion 

velocity v. Each symbol represents a different type of the final cluster configuration, as shown by snapshots (middle). The 
black lines indicate an ideal hard sphere scaling, (b) Rescaled state diagram, showing that state transitions are proportional 
to vr?I(F g^') for small self-propulsion and dipole strengths, (c) Trajectories of the individual dipoles observed in the moving 
center-of-mass frame [5^. (d) Comparison of the emerging clusters for fixed reduced dipole strength m^/(ecr^) and varied 
reduced self-propulsion vr ja for an instantaneous and an adiabatic switching protocol. 


with each other (0 = 180°), leading to a configuration of 
two clusters propagating in opposing directions. If v ex¬ 
ceeds a critical threshold, a fission into several units will 
ultimately occur. Again, we can discriminate between 
a state where the detached dipoles move in the opposite 
(90° < (j) < 180°) or in the same direction (0o < 0 < 90°) 
as the backbone of the initial cluster. As a final remark, 
the hard sphere scaling has a similar performance in the 
angle (j) as shown in Fig. HJ^c). 

Finally, the center-of-mass velocity Vc = |dRc(t)/dt| 
(see Eq. ([6|)) is shown in Fig. [5] for all emerging clus¬ 
ter conformations. While the kite cluster (yellow filled 
diamond) and the split-into-two-units state (blue open 
squares) have a vanishing center-of-mass velocity Vc, it 
only vanishes temporarily during rearrangement into the 
anti-Y-junction (green filled pentagon). 


B. Compact cluster (N = 5) 

We finally study an initial three-dimensional clus¬ 
ter with a nonvanishing dipole moment for A" = 5 
dipoles. Such a compact cluster of passive dipoles is 
metastable only for w?/{ea^) < 4.5, see Fig. [bfa). For 
rin?/{e(7^) > 4.5 this cluster will spontaneously trans¬ 
form into a ring-like structure which is the correspond¬ 
ing ground state M- If self-propulsion is turned on 
instantaneously, this cluster will either stay stable and 
move on a helical trajectory as discussed in Sec. [1111 (foi* 
m?/{ea^) < 4.5) or spontaneously transform into a ring 
(for m?/{ecr‘^) > 4.5). 

The emerging state diagram for the final cluster is pre¬ 
sented in Fig. ini^a). On top of the helical cluster motion 
and the ring-like structure, fission occurs at high self¬ 
propulsion t;. As a function of increasing the compact 
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(b) 





(c) vrea^lm^ 


FIG. 4. (a) Schematic representation of the final angles 

a and 0 and (b),(c) their dependence of the self-propulsion 
velocity v, here rescaled with the fixed dipole strength. Re¬ 
spective conformations characterized by the order parameters 
are indicated by sketches. 



0.1 1 10 100 
log(i/r) 

FIG. 5. Reduced center-of-mass velocity Vc/v of an ini- 
tial Y-junction as a function of time for fixed instantaneously 
applied self-propulsion velocity v for each emerging cluster 
conformation using the color coding of Fig.[3]and fixed dipole 
strength m^/(ecr^) = 1. 


dimensional ring. 

We have documented the fission of the cluster by mon¬ 
itoring the cluster radius R of gyration defined via 

1 ^ 

= (8) 

i=l 

Results for R for an instantaneous switching protocol are 
shown in Fig. [Tl^a). Here the inverse of R^ is plotted 
such that a cluster explosion is indicated by a vanishing 
1/R‘^. The dashed line indicates a size comparable to the 
cut-off distance Tc- A comparison of R relative to this 
line clearly reveals the transition from a compact cluster 
towards fission. For completeness we have also shown 
the cluster magnetization M in Fig.[3^b) which is a good 
order parameter for a ring-like structure as M vanishes 
there. These results confirm the full phase diagram of 
Fig. Eta). 


cluster first transforms into a ring before it splits apart. 
There is a critical drive v = Vc ^ Gajr beyond which the 
compact initial cluster looses its stability. For fixed drive 
V < Uc, an interesting reentrant scenario for increasing 
dipolar strength m occurs. There is fission at small m, 
then the ring-structure is emerging as a final state, sub¬ 
sequently the compact cluster is the final state and then 
the ring structure is getting stable again. 

The comparison of different protocols is shown in 
Fig. ini^b). In qualitative accordance with the case of 
N = 4 particles discussed above there are fewer (in this 
case only two) final states. The fission is missing com¬ 
pletely for adiabatic switching and the final cluster con¬ 
formation for large self-propulsion velocities u is a two- 


V. CONCLUSION 

We have studied the dynamical response of initially 
passive and metastable soft-sphere dipole clusters to the 
onset of an internal self-propulsion. The latter is mod¬ 
elled by an internal effective force along the dipole mo¬ 
ment and the dynamics is completely overdamped in a 
solvent at low Reynolds number. Even though the cluster 
is small, a wealth of different types of motion is obtained, 
which depend on the interaction parameters, the strength 
of the self-propulsion, and the self-propulsion switch-on 
protocol. Moving cluster structures emerge which are not 
stable in equilibrium. Interestingly, if the self-propulsion 
is very large and applied quickly (“instantaneous switch¬ 
ing” ), the cluster shows a permanent fission but there is 
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FIG. 6. (a) Emerging state diagram for an initially compact cluster of = 5 dipoles for varied reduced dipole strength 

m?/{ea^) and reduced self-propulsion velocity vt ja. As in Fig. [S] each symbol represents a type of final cluster configuration, 
shown through snapshots, (b) Comparison between instantaneous switching and adiabatic switching of the self-propulsion 
velocity v for fixed dipole strength m ?= 1. 


(a) 


1.5 


0.5 


fission I 


cluster 


10 


m?‘l{ea^) = 1 —■ 

m?/{ea^) = 2.25 — ■ 
m^/(ecr^) =4 - ■ 

m^/(ecr^) = 5.76 -> 


20 


VT ! a 


30 


(b) 


0.5 





compact 
cluster 


10 


ring-like 

cluster 



20 


30 


VT ! a 


FIG. 7. (a) Inverse squared radius of gyration 1/E^ and 

(b) reduced total dipole moment Mjm of an active cluster of 
iV = 5 dipoles for fixed dipole strength and varied self¬ 
propulsion velocity v. The dashed line in (a) shows the inverse 
squared cut-off radius of the WCA potential l/r^. 


no fission if the drive is applied slowly (“adiabatic switch¬ 
ing”). The center-of-mass of the cluster moves on linear, 
circular, or helical paths, depending on the initial mag¬ 
netization of the cluster. 

Self-propelled dipolar particles can be realized as Janus 
particles [b^, with an embedded dipole moment [bll-l^ 
where the self-propulsion [bsl, [bj is generated by self- 
diffusiophoresis. Cluster aggregation has indeed been 
observed [55|. Furthermore, dipolar particles driven by 
external magnetic fields are conceivable as experimental 
realizations at microscopic [65| as well as macroscopic 
length scales [66| . These particles can readily be observed 
in real-space such that our predictions are, in principle, 
verifiable. Another realization of self-propelled dipoles 
are active droplets which are filled with a liquid crystal. 
These possess an inner topological defect, shifted from 
their center of mass, which induces an electric dipole mo¬ 
ment. Here, the strength of the dipole moment is pro¬ 
portional to the self-propulsion, which may lead to new 
interesting effects [^, . 

Our work can and should be extended towards several 
directions for future research. First of all, soft spheres 
with a non-central dipole moment have been consid¬ 
ered recently in the passive case and the stable cluster 
structure was found to be different to that of central 
dipoles [TqI, [52|. Moreover, one can imagine that the two 
directions of dipole moment and self-propulsion are not 
collinear which is expected to lead to even more com¬ 
plex dynamics. In this case, fission is expected also for 
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slow switching. For microswimmers, Brownian fluctua¬ 
tions induced by the solvent needs to be considered and 
incorporated into the dynamics HSlsi. 

A more detailed modelling which we shall pursue in the 
future is an inclusion of hydrodynamic interactions by us¬ 
ing more complicated friction tensors. The friction tensor 
corresponding to the direct forces (or body forces) is de¬ 
rived from the pairwise interparticle potential. It can be 
treated by the Oseen or Rotne-Prager tensor [tqI-It^ , as 
would be done for passive particles. On this level, strictly 
speaking a dumbbell has another friction tensor than of 
two spheres. It is important to note that these friction 
tensors do not affect the initial equilibrium structure of 
the cluster. The swimming process, however, needs to 
be modelled with a tensor that decays faster with inter¬ 
particle distance than the Oseen tensor Izilll. We 
expect that solving the coupled equation of motion could 
give rise to new unexpected cluster dynamics as both in¬ 
teractions (the hydrodynamic and the dipolar one) are 
long-ranged and therefore compete. Moreover, the be¬ 


haviour will depend on the hydrodynamic boundary con¬ 
ditions. An unbounded solvent around the cluster will be 
described by a different friction tensor than the motion 
of dipolar particles on a substrate or on a pending 
air-liquid interface which will make the full problem even 
more complicated. Finally, particle shapes different from 
that of a sphere can be studied like C- or L-shaped parti¬ 
cles [^, [7y, UM 5 which tend to perform a circular motion 
if they are self-propelled even in absence of any dipolar 
interaction. 
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